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We report on local magnetization, tunnel diode oscillator, and specific heat measurements in a 
series of Ba(Ni I Fei_ a; )2As2 single crystals (0.26 < x < 0.74). We show that the London penetration 
depth A(T) = A(0) + AA(T) scales as A(0) oc 1/T C °' 85±0 - 2 , AA(T) oc T 2 3±0 3 (for T < T c /3) and 
dAX/dT 2 oc l/T 2 ' 8±0 ' 3 in both underdoped and overdoped samples. Moreover, the slope of the 
upper critical field (H' c2 = —(dH C 2 / dT)\ T ^> Tc ) decreases with T c in overdoped samples but increases 
with decreasing T c in underdoped samples. The remarkable variation of A(0) with T c and the non 
exponential temperature dependence of AA clearly indicates that pair breaking effects are important 
in this system. We show that the observed scalings strongly suggest that those pair breaking effects 
could be associated with quantum fluctuations near 3D superconducting critical points. 



I. INTRODUCTION 

The discovery of superconductivity up to 55K in iron- 
based pnictides [I] has generated tremendous interest. 
Even though the symmetry of the order parameter in 
this multi-band system is not yet determined with cer- 
tainty, a popular model is based on magnetic fluctua- 
tions associated with a sign reversal of the order param- 
eter between the hole and electron sheets of the Fermi 
surface (so called s± model 0). It has then been sug- 
gested by V.G.Kogan [31 0] that the critical temperature 
(T c ) could be strongly suppressed not only by scattering 
breaking the time reversal symmetry (spin-flip scatter- 
ing) but by any scattering mechanism. As a consequence, 
for T c « T c .o (Tcfi being the critical temperature in the 
absence of scattering) and an average of the order pa- 
rameter over the Fermi surface being equal to zero, the 
superfluid density (p s ) is expected to vary as: 

Ps oc 1/A 2 oc T 2 - T 2 (1) 

leading to a London penetration depth (A) scaling as : 
A(0) oc 1/T C and AA(T) = A(T) - A(0) oc T 2 /T 3 (for 
T — > 0). Moreover, the slope of the upper critical field 
close to T c , H' c2 = —(dH c2 /dT)\ T ^ Tc is then expected to 
be proportional to T c and finally the specific heat jump 
at T c , AC P oc T 3 . Some indications for the AA(T) [5] or 
H' c2 and AC P [6] scalings have been observed in various 
pnictides but a systematic analysis of the evolution of all 
those quantities on a given system was still lacking. 

Moreover, the dependence of p s (0) on T c remains 
controversial. It has been initially suggested [7] that 
p s (0) could scale as T c but, measurements in overdoped 



Ba(Co a: Fei_ a; )2As2 rather suggested that p s (0) could be 
either proportional to X" 2 8J or, on the contrary, almost 
T c independent [9. . Those later measurements also in- 
dicated a strong reduction of the superfluid density in 
underdoped samples, which has been attributed to the 
coexistence of superconductivity and magnetism. 

It is worth noting that somehow contradictory re- 
sults have also been observed in cuprates, emphasizing 
the fact that the p s (0) vs T c dependence can be very 
sensitive to the sample quality and/or dimensionality 
[TO] . As in pnictides, initial measurements suggested 
that p s (0) oc T c in samples close to optimal doping [TTJ 
and this dependence found a straightforward explana- 
tion in this quasi 2D system assuming that T c is close 
to the Kosterlitz-Thouless-Berezinskii transition temper- 
ature T c ~ Tktd = ^ot/8np ^ 2 (t being the interlayer 
spacing). However, strong deviations from this behav- 
ior have been observed in highly underdoped samples for 
which p s (0) was found to scale as T£ with 1.6 <( < 2.3 
in thick films [H] and ultraclean crystals [131 [2] or £ <~ 1 
in very thin films [TU] . This change in £ from ~ 2 to 1 has 
been attributed to a dimensional crossover associated to 
the proximity of a quantum critical point [15] . 

We present here a detailed analysis of the doping de- 
pendence of T c , A(0), AA(T) and specific heat in the 
Ba(Ni x Fei_ x )2As2 system. The measurements have been 
carried out in a series of Ba(Ni E Fei_ x )2As2 single crys- 
tals grown by self flux method. Details on the sam- 
ple elaboration are given in [16]. We will show that 
A(0) oc l/T c °- 85±0 - 2 , AA(T) a T 2 3±0 - 3 (for T < T c /3) 
and dAX/dT^ Q oc l/T 2 - 8±0 - 3 in both underdoped and 
overdoped samples. Moreover, H' c2 decreases with T c in 
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overdoped samples but increases with decreasing T c for 
underdoped samples. Those results strongly suggest the 
presence of pair breaking effects and we will show that 
they can be consistently described assuming that those 
pair breaking effects are associated with the proximity of 
superconducting quantum critical points. 



II. LONDON PENETRATION DEPTH 

A. Lower critical field measurements 

The local field has been measured by placing the sam- 
ples on miniature GaAs-based quantum well Hall sensors. 
The external field was increased up to H a and swept back 
to zero in order to measure the remanent field (B rem ) 
trapped in the sample after the field cycle. All measure- 
ments were performed for ff a ||c. In the Meissner state, 
the external field is fully screened out and B rem is equal 
to zero. Vortices start to penetrate into the sample for 
H a > Hp (the first penetration field) and remain par- 
tially pinned in the sample when the field is swept back 
to zero, leading to a finite B rem value (see Fig. la). For 
H a > Hp, the remanent fields increases approximatively 

2.5 and H p as hence 
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as (H a - Hp) 1 '-* with l/ 7 - 1.5 
been obtained by linearly extrapolating B^ em to zero (see 
Fig. la). In order to avoid spurious effects associated to 
strong pinning preventing the vortex diffusion to the cen- 
ter of the sample (and hence the absence of any signal 
on probes located close to the center of the sample) H p 
has been measured on several locations with an array 
of miniature probes. As expected the measured H p value 
(slightly) increases as the distance between the probe and 
the sample edge increases. The H p values reported in Ta- 
ble 1 have been obtained with probes located at ~ 5 to 
10/im from the edge and we estimate to ~ 20 % the pos- 
sible overestimation of this field. 

This first penetration field is smaller than the lower 
critical field (H c \) due to the expulsion of the flux lines 
for H a < Hp which increases the local field in the 
vicinity of the sample edges so that H c i — aH p with 
a > 1. In the presence of geometrical barriers (GB) 
[FT] , cxq b ~ tanh^d/w) -1 / 2 where /3 varies from 0.36 in 
strips to 0.67 in disks (d and w being the thickness and 
width of the sample, respectively) whereas for elliptical 
samples (i.e. without geometrical barriers) the standard 
"demagnetization" factor a e ui P tical — 1/(1 — N) cx w/d. 
All samples were chosen to present very similar aspect 
ratios (see Table 1) leading to very similar corrections 
for the whole series whatever the origin of a. In the 
following, we approximated the samples by disks lead- 
ing to 2.4 < ctGB < 2.8 from one sample to the other 
(see Table 1). However, it is important to note that 
a standard "demagnetization" correction would lead to 
2.7 < 1/(1 — A) < 3.3 for our d/w values (~ 1/4) and 
would hence only lead to an overestimation of H c \ by 
~ 12% i.e. an underestimation of A of ~ 6%. In the fol- 
lowing, to evaluate the error bars, the uncertainty on a 
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FIG. 1. (a) i?7em as a function of the applied field for the 
indicated temperature in Ba(NL c Fei_ a; )2As2 with x = 3.4% 
(B rem being the remanent field trapped in the sample after 
a field excursion up to H a - see text for details - and 7 is an 
exponent ~ 0.4 — 0.6). The first penetration field is deduced 
from the linear extrapolation of B^ em to zero. Inset (color on- 
line) : Same as main panel for different x values at T — 1.3K 
(the different curves are shifted vertically by 0.2 for clarity), 
(b) Temperature dependence of the lower critical field H c \ 
in Ba(Nia;Fei_^)2As2 single crystals (for the indicated dop- 
ing concentrations x). The solid lines are the values expected 
from Eq.(3) introducing only one free parameter (correspond- 
ing to the absolute value at zero T) for the whole set of data 
(see text for details). The shaded cones correspond to the 
values expected from the amplitude of the specific heat jump 
at T c (see Fig. 4) and the slope of the upper critical field (see 
Fig. 5 and text for details). 



is taken to ±20%. The temperature dependences of H c \ 
(for the indicated x values) are displayed in Fig. lb. 

Finally, A(0) has been deduced from the upper and 
lower critical fields writing : [IqH c2 = <&q/2-k£ > 2 (with 
H c2 {0) ~ 0.7.T c .(dH c2 /dT) lT ^ Tc [IE], see section III) 
and \iqH c \ — ($ /47rA 2 ).(Ln(K) + c(k)) where k = A/£ 
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TABLE I. x, d, w and I are the composition, thickness, width and length of the Ba(Ni a: Fei_ a ;)2As2 samples, respectively. 
The first penetration field H p has been defined as the field above which a remanent field is observed in the sample and the 
lower critical field H c \ (see Fig.l) is deduced from H p introducing demagnetization corrections due to geometrical barriers 
(noH c i[G] — acB x fj,oH p [G] (see text). A is the London penetration depth (the values between brackets are deduced from 
specific heat measurements) and /3 = dX/dT 2 (see Fig. 4). H C 2 is the upper critical field (see Fig. 5, the values deduced from 
the irreversibility field are marked by an *). All values are given for T — > 0. AC P /T C is the amplitude of the specific heat 
jump at T = T c (in mJ.moPVK -2 ). G; is the Ginzburg number (see text) and T c the critical temperature. Measurements 
techniques : HP=Hall Probe, Cp=Specific Heat and TDO=Tunnel Diode Oscillator. 
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and c(k) is a n dependent function tending towards ~ 0.5 
for large k values, k has hence been deduced from the 
H C 2/H c i ratio (see table 1) and the corresponding A(0) = 
k£(0) values have been reported on Fig. 2a as a function 
of T c . As shown, similar values (A(0) - l/T c °- 85±a2 ) 
are obtained on both sides on the superconducting dome 
in striking contrast with the result obtained by Gordon 
et al. in Ba(Ni x Fei_2;)2As2 |9J. Our measurements do 
hence not support the scenario of a sharp increase of A(0) 
due to the presence of a magnetic gap in underdoped sam- 
ples but rather suggest that only a small fraction of the 
Fermi surface is affected by the antiferromagnetic cou- 
pling. 

Note that (dH cl /dT)\ T ^ Tc (and hence A(0)) can also 
be deduced from specific heat measurements (see section 
III). Indeed, the amplitude of the specific heat jump 
at T c is equal to : AC P = dj, T c ).(dH c /dT)? T ^ Tc ~ 
(fi T c ).(dH c2 /dT) lT ^ Te .(dH cl /dT) lT ^ T J[ln(K + 6.5]. 
Deducing H C 2(T) from the shift of the anomaly under 
magnetic field (see section III. A), one directly obtains 
(dH c i/dT)\T^T c an d X Cp (0) can then be calculated 



using the same procedure as above, assuming that 
H cl (0) - 0.7 .T c .(dH cl / 'dT) ]T ^ Tc . The corresponding 
H cl (T -> T c ) = (T — T c ) x (dH c i/dT) ]T ^ Tc and A c *(0) 
values have been reported on Fig. lb (shaded cones) and 
Fig. 2a (crosses), respectively. It is important to note 
that those values do not depend on any demagnetization 
correction. As shown in Fig. 2a, the as-deduced X Cp 
values are in reasonable agreement with those deduced 
from Hp (X Hp ) measurements (X Hp /A Cp ~ 0.7) and both 
measurements lead to very similar dependences on T c . 



The most remarkable feature is the very large variation 
of the superfluid density which varies over one order of 
magnitude for T c values ranging from ~ 6.8K to ~ 20. 2K. 
Such a strong dependence is unexpected in conventional 
superconductors, but strongly suggests the presence of 
pair breaking effects. 
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FIG. 2. (a) London penetration depth at T — >■ (A(0)) as a 
function of the critical temperature T c in Ba(Ni I Fei_ :!; )2As2 
crystals deduced from H c i measurements (squares, see Fig.l) 
and specific heat measurements (crosses, see text for details). 
The thick line is a guide to the eyes varying as 1/TJ? .85. (b) 
Slope of the AA vs T 2 curve (see Fig. 3) as a function of T c 
in samples from the same batch deduced from TDO measure- 
ments (squares). The crosses correspond to X(0)/2T 2 (see 
text for details) introducing the A(0) values deduced from C p 
measurements (see upper panel). 



B. Tunnel Diode Oscillator Technique 

The samples were glued at the end of a sapphire rod 
which was introduced in a coil of inductance L. The 
variation of the London magnetic penetration depth in- 
duces a change in L and hence a shift of the resonant fre- 
quency Sf(T) = f(T) — f(T m i n ) of a LC oscillating circuit 
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FIG. 3. London penetration depth X mix = \ a b + (d/w)X c (~ 
Aai, for w >> d) as a function of (T/T c ) 2 in Ba(Ni :E Fei_ :E )2As2 
single crystals for the indicated x values. As in Fig. 2, closed, 
half-closed and open symbols refer to underdoped, optimally 
doped and overdoped samples, respectively. Inset : best fit 
values for the exponent n, assuming that A A oc T n . 



(14MHz) driven by a Tunnel Diode. This shift, renor- 
malized to the one corresponding to the extraction of the 
sample from the coil (A/o) is then equal to the fraction 
(SV/V) of the sample which is penetrated by the field. 
For H\\c, 5V is related to the in-plane penetration depth 
X a b through some calibration constant depending on the 
sample geometry. However, this constant can be altered 
by edge roughness effects (see discussion in [3DJ) and we 
have hence decided to perform all measurements with 
H\\ab. Indeed, the surfaces parallel to the a6-planes are 
much flatter and 5V/V is, in this case, directly given by 
SV/V ~ 2(A c /u>+A ab /d) - 5V/V ~ 2/dx [\ ab + (d/w)\ c ] 
without any geometrical correction (A c being the pen- 
etration depth parallel to the c-axis). In contrast to 
H c i measurements for which we used rather "thick" 
(d/w ~ 1/4) samples in order to reduce the uncertainty 
related to geometrical corrections, we have, in this case, 
selected very thin samples (d/w << 1) (see Table 1) so 
that A a £> + [d/w)\ c = X m i x ~ X a t (for weakly anisotropic 
systems [!§]). 

The temperature dependence of the penetration depth 
is clearly non exponential in all measured samples but 
can be well described by a power law : AX mix (T) — 
X m ix(T) - X mix (0) = AT~ 2 - 3±0 - 3 for T < T c /3 (see 
Fig. 3). A very similar behavior has been reported in 
a large number of pnictides in both 1111 and 122 sys- 
tems [S] (as well as in Fe(Se,Te), see [2D] and references 
therein). This dependence is another clear indication for 
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FIG. 4. Temperature dependence of the field dependent 
part of the specific heat for the indicated field values in 
Ba(Ni x Fei_ x ) 2 As 2 . 



the presence of pair breaking effects but it is also im- 
portant to note that the slope of the AA vs T 2 curve, 
dAX/dT 2 T ^ Q is proportional to 1/T C 2 8 (see discussion in 
section IV. A). 
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FIG. 5. Temperature dependence of the upper critical field 
H C 2 deduced from the mid-point of the specific heat anomaly 
(squares, see also Fig. 2 for symbol correspondance). The 
crosses correspond to the irreversibility field (Hi rr ) deduced 
from the onset of diamagnetic screening and circles to H C 2 
values deduced from Hi rr (see text for details) . 



titative C p values, special care has been taken in the cali- 
bration procedure to measure a copper standard and the 
addenda in the exact same conditions. For x = 3.4, 4.2 
and 6.5 % well defined specific anomalies were obtained 
in zero field (see Fig. 4) and this anomaly progressively 
shifted with magnetic field. T c2 has hence been defined 
as the temperature corresponding to the mid point of 
the transition for a given value of the external field. The 
transition became broader and weaker for x = 2.9 % (not 
shown) and finally too weak (and/or broad) to allow any 
accurate determination of H c2 from specific heat mea- 
surements for x = 2.6 % and x = 7.4 %. The H c2 {T) 
values have been reported in Fig. 5 and the correspond- 
ing H' c2 values are displayed in Fig. 6b together with the 
T c versus x data. 



III. UPPER CRITICAL FIELD 

A. Specific heat measurements 

Finally, C p measurements have been performed in 
magnetic fields up to 28 T using an high sensitivity AC 
technique (typically 1 part in 10 3 ). Heat was supplied 
to the sample by a light emitting diode via an optical 
fiber and the corresponding temperature oscillations were 
recorded with a thermocouple. In order to obtain quan- 



B. Transmittivity measurements 

In order to obtain the upper critical field for the two 
low T c samples, we performed transmittivity measure- 
ments (T' H ). The AC component of the local induction 
on the Hall probe was recorded in presence of an AC mod- 
ulation field (h ac ~ 1G, u! ~ 200 Hz) and T' H has been 
defined as '. Tu — t-tj trS^ 2-* \ o ~ ^ ^rl w . The irre- 

H (Bac{T>T c )~Bac(T«T c )) 

versibility field (H irr ) was defined as the onset of diamag- 
netic response (i.e. locus of minimal observable screen- 
ing current). Writing the Ginzburg-Landau free energy 
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functional in terms of Lowest Landau Level (LLL) eigen- 
functions, H C 2 can then be deduced from Hi rr through 
ED : 

(1 - h)(l - t 2 ) 1/3 = (th) 2/3 x CGi 1/3 (2) 

where t = T/T c , h = H trr (T)/ H c2 (T), C is a constant 
depending on the amount of disorder present in the sam- 
ple and Gi = | [fcsT c /£oe£(0))] 2 is the Ginzburg number 
(see Table 1) with e = ($ /47rA(0)) 2 the vortex line 
energy and e the anisotropy [19] . 

As shown in Fig. 5, the H C 2{T) values (diamonds) de- 

1 /3 

duced from H irr (T) (crosses) taking CG i ~ 1 are in 
very good agreement with the values deduced from the 
specific heat measurements (squares) for both x = 3.4 
and x = 6.5 %. Note that the irreversibility line lies sig- 
nificantly below the H C 2(T) line, clearly suggesting the 
presence of vortex liquid phase in this systems |22| . The 
corresponding C value (~ 15) is very close to the one 
obtained in (Ba,K)Fe2As2 [23] and we have calculated 
the H c2 values for x = 2.6 and 7.4 % using this C value 
(i.e. taking CG 1 / 3 ~ 1.1 ± 0.2 and CG 1 / 3 ~ 0.8 ± 0.2 for 
x = 2.6 and x = 7.4, respectively) We hence confirmed 
the decrease of H' c2 with T c in overdoped samples and 
the increase of H' c2 for T c — > for low doping contents 
(see also Fig. 6b). Note that a similar behavior has been 
reported by Vinod al. [24] in Ba(Ni 2; Fei_ a ;)2As2. Even 
though the determination of H C 2 from transport measure- 
ments may be altered by fluctuation effects, those mea- 
surements also suggested that H' c2 could be proportional 
to T c in overdoped samples and that the H' c2 /T c ratio in- 
creases in underdoped samples. Finally, note also that a 
different scaling has been recently obtained in overdoped 
BaFe2(As,P)2 samples in which H' c2 rather scales as T 2 
[25] (see discussion in section IV. B). 



IV. DISCUSSION AND CONCLUDING 
REMARKS 

In summary, we have shown that : 

i) the penetration depth A strongly varies with the crit- 
ical temperature of the sample, 

ii) A(0) scales as Tjp- 85±0 - 2 for both underdoped and 
overdoped samples, 

iii) the temperature dependence of A is non exponential 
but varies as AA oc T 23±03 for T < T c /3 , 

iv) the slope of the upper critical field decreases with 
decreasing T c in overdoped samples but increases with 
decreasing T c in underdoped samples. 

Points i) and iii) strongly suggest that pair breaking ef- 
fects are important. As pointed out by V.G.Kogan [2l|4], 
in the s± model, the Cooper pairs are expected to be 
be very sensitive to all scattering events and the super- 
conducting condensate is hence progressively destroyed 
on both side of the superconducting dome, in agreement 
with the observation of a non residual Sommerfeld co- 
efficient in specific heat data in Co-doped samples [50] . 
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FIG. 6. (color online) (a) 1/£A (closed symbols), T c 15 (small 
open symbols) and T C (AC* P /T C ) - 5 (large open symbols) as a 
function of doping in both Ba(Fe a: ,Nii_ a ;)2As2 (present work, 
(black) circles) and Ba(Fe :E ,Coi_ 2 ,)2As2(blue) squares). The 
x values in Co-doped samples have here been rescaled by 
a factor ~ 0.6. (b) Left scale : slope of the critical field 
H' c2 — dH C 2/ 'dT\ T ^ T< , (solid lines in Fig. 5b) as a function 
of the composition (x) in Ba(Ni I Fei_ 2 :)2As2 single crystals 
where H C 2 has been deduced from specific heat measure- 
ments (open squares) and from the irreversibility field (solid 
squares). Right scale : critical temperature T c (solid circles) 
as a function of x. 



However, it is still unclear whether the system lies close 
to the critical disorder for which T c — > and whether 
the gap vanishes or not. Alternatively, those pair break- 
ing effects can be associated with quantum fluctuations 
[IB] . Those fluctuations require a finite value of the cou- 
pling strength even for T c — > [28]. They can hence 
not be observed in a standard BCS superconductors but 
can in the presence of pair breaking effects as supercon- 
ductivity only develops above some finite value of the 
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coupling strength in this case. The scaling properties 
associated with those two scenarios (strong pair break- 
ing effects close to critical disorder and superconducting 
quantum fluctuations) are discussed below. 

A. Strong pair breaking effects for < SI >= 

As discussed by V.G.Kogan [3JH], the supernuid den- 
sity is expected to vary as (T c (x) 2 — T 2 ) (Eq.(l)) for an 
average of the order parameter over the Fermi surface 
< Q >= and vanishingly small critical temperatures 
(see also [H]). Writting 1/A(T) 2 = 1/(A(0) + AA(T)) 2 ~ 
(1/A(0) 2 )x(l-2AA(T)/A(0)) (for AA << A(0)) one then 
expects that AA(T) - f3T 2 with /3 = A(0)/2T 2 . The ex- 
ponent n slightly larger than 2 observed experimentally 
can be attributed to a non zero gap value but in order 
to compare the different samples within each other, we 
assumed in the discussion below that n — 2. Neglecting 
the - small - variation of ln(n) with x one also expects 
that : 

H cl (x,T) = K{x) x (T c (x) 2 -T 2 ). (3) 

This equation leads to a very satisfying agreement to the 
data introducing only one single adjustable parameter 
(A^a;) = K ) for the whole set of x values (see solid 
lines in Fig. lb). As mentioned above, the aim is here not 
to obtain the best fit to the data (the agreement could 
for instance be improved by adjusting the exponent n) 
but to emphasize that the whole set of data can be re- 
produced in a very satisfactory way only assuming that 
Eq.(l) is valid. Eq.(l) also suggests that A(0) cx 1/T C 
in very reasonable agreement with the A(0) = A/T®- 85 
scaling observed in Fig. 2a (point ii)). One then expects 
P = A(0)/2T 2 = B/T 2m in very good agreement with 
the experimental data (see Fig. 2b). TDO measurements 
hence also lead to A(0) ~ A/T - 85 . The corresponding A 
value is on the order of ~ 5600 nmK - 85 e.g. about two 
times larger than the one obtained from H c i measure- 
ments (A ~ 2800 nmK - 85 ) but in very fair agreement 
with the one deduced from C p measurements (~ 4000 
nmK 85 , see corresponding A(0)/2T 2 values in Fig. 2b). 
It is important to note that, this scaling law has been 
obtained through three independent techniques (with an 
average A value ~ 4100 nmK 85 ). Note that much larger 
(3 values (from ~ 1.2 nm/K 2 in the optimally doped sam- 
ple to ~ 45 nm/K 2 for x = 7.2%) have been obtained by 
Martin et at 29J in Ni doped crystals but those val- 
ues would correspond to A(0) ranging from ~ 900 nm 
to ~ 1600 nm (in the pair breaking scenario) i.e. much 
larger than any experimental value. 

The scaling properties of the magnetic penetration 
depth are hence in very reasonable agreement with the 
strong pair breaking model suggested by V.G.Kogan 
[3j 0] and it is reasonable to attribute, in this scenario, 
the small deviations from the expected scaling laws to 
the non zero gap value. However, it is important to note 
that this model assumes a very strong reduction of T c 



suggesting that T Cj0 could largely exceed 100 K in pnic- 
tides in the absence of scattering which remains very puz- 
zling. Moreover, sample dependent parameter (such as 
the density of states or the Fermi velocity) enter in the 
prefactors of the scaling functions and it is hence quite 
surprising to obtain one single A value (i.e. K value for 
the H c i data) for all doping contents [37] . Moreover, even 
though V.G.Kogan pointed out in [1] that the H' c2 cx T c 
dependence expected for < ft >= would actually be re- 
versed (H' c2 cx 1/T C ) for < Q >^ 0, the strong difference 
in H' c2 (T c ) between underdoped and overdoped samples 
(point iv)) and not in A(0) vs T c , can hardly be explained 
in the framework of this model. 



B. Proximity of a superconducting quantum 
critical point 

In the presence of quantum fluctuations, static and 
dynamical properties are inextricably mixed, so that 
the value of the dynamic exponent z directly enters in 
the scaling properties of the thermodynamical properties 
[15] . However, for d=2 (d being the dimension), z cancels 
out in the p s (0) vs T c scaling and p s (0) cx T c for all z. 
On the other hand, for d=3, d + z > 4 (i.e. larger than 
the Gaussian end point value) and the superfluid density 
is expected to vary as : p s (0) cx \p — p c \ where p is the pa- 
rameter driving the transition (here the doping content) 
and p c its critical value. The critical temperature varies 
as : T c cx \p — p c |* where the shift exponent \& is related 
to z through ^> = z/(z + 1) 45j and one expects : 

p.(o) cx rj* +1 )/* (4) 

with z = 2 in standard superconducting quantum critical 
points and z — 1 in the presence of nodes in the supercon- 
ducting gap [28]. One then expects either p s (0) cx T 15 
for z — 2 or p s (Q) cx T 2 for z — 1. The change from 
p s (0) cx T c to p s (0) cx T 2 with the sample thickness in 
underdoped YBaCuO has hence been interpreted as a di- 
mensional crossover from 2D thin samples to 3D samples 
with 2 = 1 [TP] . 

In our case, p s cx T^ 7±0A and it is hence impossi- 
ble to distinguish between the z = 1 and z = 2 cases 
(d = 3). On the other hand, 1/£q is also expected to scale 
as Tc and, as for the strong pair breaking model dis- 
cussed above, it would be difficult to explain the different 
H' c2 (T c ) dependences observed for underdoped and over- 
doped samples assuming that they are in the clean limit 
(H' c2 cx 1/£qT c ). This difference can however be under- 
stood in the dirty limit for which £ ~ (£o0°' 5 where I and 
fo are the mean free path and coherence length without 
disorder, respectively. A non symmetric l[x) dependence 
would lead to different H' c2 dependence for underdoped 
and overdoped samples. Note that some difference would 
also be expected for A ~ Ao(£o/0°' 5 (^o being the pen- 
etration depth without disorder) but this difference is 
much weaker than for H' c2 and remains within our error 
bars. 



Interestingly, I is expected to cancel out in the £A prod- 
uct which is equal to £oAo f° r both clean (I > £0) an d 
dirty (I < £ ) samples and one expects : 

l/eAocTf+W* oc |z-a; c | (5) 

on both sides of the superconducting dome, indepen- 
dently of the sample quality. Note that, two quantum 
critical points are present in this scenario, correspond- 
ing to the two end points of the dome. As shown in 
Fig. 6a, 1/£A (closed circles) and T^ 2 (small open cir- 
cles) present very similar dependence on the doping con- 
tent clearly suggesting that 1/£A oc T c 3 ^ 2 i.e. that z = 2. 
Moreover, they both decrease roughly linearly with x on 
either side of the superconducting dome in good agree- 
ment with Eq.2. 

Moreover, as pointed out in section II. A, the jump 
of the specific heat at T c : AC P /T C oc [1/£AT C ] 2 and 
we have also reported on Fig.6a the T c -J AC P /T C val- 



ues (large open circles) deduced from our specific heat 
measurements. As shown, this quantity consistently fol- 
lows the expected behavior. Moreover, one can similarly 
plot 1/fA, T c 15 and T cy / AC P /T C as a function of doping 
for Ba(Fe 2 Coi_ 2; )2As2 samples (squares) taking A values 
from [5] and [31], H' c2 data from (24] and T cy fAC p /T c 
values from [5J. As shown, very similar scalings are ob- 
tained in both systems (the x values in Co-doped samples 
have here been rescaled by a factor ~ 0.6 to take into ac- 
count the difference in the electronic valency between Co 
and Ni atoms) emphasizing the strong similarity between 
the two systems. 
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